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ABSTRACT 

Interacting galaxies are well-known for their high star formation rates and rich star cluster 
populations, but it is also recognized that the rapidly changing tidal field can efficiently de- 
stroy clusters. We use numerical simulations of merging disc galaxies to investigate which 
' mechanism dominates. The simulations include a model for the formation and evolution of 

the entire star cluster population, accounting for the evaporation of clusters due to two-body 
relaxation and tidal shocks. We find that the dynamical heating of stellar clusters by tidal 
shocks is about an order of magnitude higher in interacting galaxies than in isolated galaxies. 
This is driven by the increased gas density, and is sufficient to destroy star clusters at a higher 
rate than new clusters are formed: the total number of stellar clusters in the merger remnant 
| is 2-50% of the amount in the progenitor discs, with low-mass clusters being disrupted pref- 

erentially. By adopting observationally motivated selection criteria, we find that the observed 
surplus of star clusters in nearby merging galaxies with respect to isolated systems is caused 
by the observational bias to detect young, massive clusters, and marks a transient phase in 
£N) ■ galaxy evolution. We provide a general expression for the survival fraction of clusters, which 

increases with the gas depletion time-scale, reflecting that both the formation and the destruc- 
tion of clusters are driven by the growth of the gas density. Due to the preferential disruption 
of low-mass clusters, the mass distribution of the surviving star clusters in a merger remnant 
develops a peak at a mass of about 10 3 M Q , which evolves to higher masses at a rate of 0.3- 
0.4 dex per Gyr. Briefly after a merger, the peak mass depends weakly on the galactocentric 
/\ ' radius, but this correlation disappears as the system ages due to the destruction of clusters on 

\ eccentric orbits. We discuss the similarities between the cluster populations of the simulated 

merger remnants and (young) globular cluster systems. Our results suggest that the combi- 
nation of cluster formation and destruction should be widespread in the dense star-forming 
environments at high redshifts, which could provide a natural origin to present-day globular 
cluster systems. 

Key words: galaxies: evolution - galaxies: interactions - galaxies: star clusters - galaxies: 
starburst - galaxies: kinematics and dynamics - (Galaxy:) globular clusters: general 



1 INTRODUCTION 



Merging and interacting galaxies host huge star bursts and large 
populations of young massive stellar clusters (e.g. iHoltzman et al.l 
1 19921 : ISchweizeret~allll996l : IWhitmore et"ai]|l999h . A galaxy in- 
teraction triggers inflows of interstellar gas to wards the galaxv 



centres, where it fuels a burst of star formation (Hernquis 



iMihos & Hemquistl 1 19961 : H arnes & Her nquist 



1996). Merger- 



induced starbursts play a central role in the history of the uni- 
verse, as galaxies ar e thought to have formed through hierar- 
chical merging ( e.g. IWhite &Reesl 1 19781 : IWhite & Fren3 1 199 it 
ICole et ail |2000). Some fracti on of this star formation takes 
place in compact ste llar clusters telmegreenlll983l ;l Whitmore et al.l 
ll999l:lBressert et alj|20ich w ith masses in the range 10 2 -10 8 M 
jPortegies Zwart et alj 120100 . The clusters that remain after a 
merger are often used as foss ils to trace the formation history of 
the galaxy jLarsen et alj200lT) . 
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During the past two decades, observations with the Hubble 
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Space Telescope have revealed that many nearby ongoing galaxy 
mergers host exceptionally rich star cluster popul ations with clus- 
ter masses exceeding 10 7 Mm iSchweizerlfl982l: iHoltzman et alJ 
19921 : lMilleretal.ll 19971 ; ISchweizer & Seitzerll 19981 ; iBastian etafl 



_to the pe 

lar medium (ISM) JSchweizeg[l987l ; Ushman & Zepj|l992l) . The 
multitude of star clusters suggests that they are useful tracers of 
past galaxy mergers, especially because they are easily observed 
up to distances of several tens of megaparsecs. The observed clus- 
ters (> 10 4 Mg) are distributed acc ording to a power la w with 
index —2 down to the detection limit (Zha ng &Falllll999h . These 
clusters are thought to be just the 'tip of the iceberg', since the ini- 
tial cluster mass function (ICMF) appears to continue beyond the 
det ection limit and down to a ce rtain physical lower mass limit (see 
e.g. lPortegies Zwart et all2010l) . 

However, high gas densities and tidal shocks, both of which 
are prevalent in coalescing galaxies, are known to have a disruptive 
effect on star clusters dSpitzen[l958l ; IWeinberd [19941 ; iGieles et all 
2006). The destruction rate of star clusters decreases with increas- 
ing cluster mass and density ( Spitzeil l 19871 : lLamers et al] |2005 a ) Q 
This indicates that the effects of star cluster disruption could be 
masked by observational selection effects and go unnoticed in ob- 
servations, i.e. the brightest and therefore most massive clusters are 
easiest to detect but also least affected by disruption. 

If the ICMF is universal, i.e. all stellar clusters are formed ac- 
cordin g to a power law with ind ex —2 throughout space and time 
(e.g. iKruiissen & Cooped 1201 lh . then the important role of clus- 
ter disruption is supported by the old ('globular') star cluster sys- 
tems that are observed in nearby spiral and giant elliptical galaxies, 
which are strongly lacking low-mass clusters wi th respect to the 
young populations i n presently mergi ng galaxies 1 Vesperinil l200ll ; 
iFall & Zhan^l200ll ; lElmegreerJl20ld) . For a power law ICMF, the 
size-of-sample effect would also require that the most massive clus- 
ters are formed in the largest bursts of star formation, implying that 
globular clusters originate from starburst environments. The ques- 
tion thus arises whether or not the disruption of star clusters dom- 
inates over their formation in starburst galaxies. This is not easily 
determined on analytical grounds. 

Globular cl uster systems are present over most of the galaxy 
mass range (e.g. IPeng et alj [2008). and as such it is evident that 
they were not only formed in interactions between massive spiral 
galaxies; the presence of globular clusters in dwarf galaxies sug- 
gests that these also endured starbursts during their early evolution. 
While the globular clusters of dwarf galaxies are generally metal- 
poor, the colour distribution of globular clust ers is often bimodal in 
massive spiral galaxies and giant ellipticals JSearle & Zin 



iForbes et al.l 1 19971; iKundu & Whitmord 1200 ll ; IPeng et alJ l2006h . 



1978; 



This colour bimodality may translate into a metallicity bimodal- 
ity, although it has recently been suggested that it i s a relic of 
a non-linear relation between colour and metalli city ( lYoon et ail 
120061 : IChies-Santos et all 1201 ll; lYoon et alJ 1201 lh . Regardless of 
whether the metallicity distribution is bi modal, a popular explana- 
tion for the broad range in metallicities dMuratov & Gnedinl fco 1 ) 
is that the metal-poo r clusters preferential ly originate from ac- 
creted dwarf galaxies (Priet o & Gnedinl2008l) . while the metal-rich 



1 Unless the environment in which they reside is so disruptive that 
it can efficiently destroy a cluster regardless of its mass. In that 
case, the rec ently argued scenario in which cluster disruption is mass- 
independent | Whitmore et alj|2007f) can arise lElmegreen & HunteJpjijTol ; 
Kruii ssenetalj|2011bl) . This would then not be universal, but depends on 
the environmental conditions. 



population was m ainly formed in-situ, eit her by disc instabilitie s 
dShapiro et ai1l20ld) or in galaxy mergers dAshman&Zepjl 19921) . 

The formation of globular clusters h as been investigated in 
several theoretical and nume ri cal studies ([ Harris & Pudritz 1994 ; 
Elmegreen & Efremovl 1 19971 ; iBekki et all 120021 : lli et aiTl2004 ; 



Bournaud et al. 2008). As expected from power law statistics, these 



studies all point to dense, gas-rich environments, which are typi- 
cally correlated with high star formation rate densities. However, 
the present-day population of globular clusters is not recovered in 
these studies, because they only concern cluster formation and con- 
tain either no description for the further evolution of clusters or a 
very simplified one. Separate studies, both analytical and numeri- 
cal, have shown that the evol ution of (globular) clusters is related 
to the galactic environment jSpitzej 19871; Baumgardt & M akino 
| 2003| : lLamers et alj2005bl:lGieles et alj2006l : Elme green & H unter 
hold : IKruiissen et alj|201 lb). To obtain a more complete, quanti- 
tative understanding of the origin of present-day globular clusters, 
it is necessary to consider their formation and further evolution si- 
multaneously. 

At present, the most commonly used method to model 
the ev olution of star cl u sters i s through JV-body s imula 



tions |Vesperini & Heggie 19971 : 



i Baumgardt & Makinol 120031: 



Praagman et al.ll2.Q10l : iRenaud et al 



Portegies Zwart et al. 
Gieles & Baumgardtl 



1998 



2008 



2011). However, this method 



is computationally too expensive to foll ow the format i on and 
evolution of the entire cluster population. IKruiissen et a D d2011bh 
therefore introduced a method in which numerical simulations 
of galaxies are supplemented with a semi-analytic model for 
the formation and evolution of star clusters, of which the results 
are consistent with (observed and simulated) formation and 
destruction rates from the literature. This model enables us to track 
the formation and evolution of the entire star cluster population 
throughout the assembly histories of galaxies. 

As a first effort to understand the (im)balance between the for- 
mation and destruction of star clusters in starbu rst environments, 
we use the method from lKruiissen et ah Il l2011bh to model the star 
cluster populations of galaxy mergers. This allows us to quantify 
the net effect of a galaxy merger on its cluster population. With this 
setup, we aim to investigate: 

(1) the relative importance of cluster formation and destruction in 
interacting galaxies; 

(2) whether galaxy mergers can produce the progenitors of 
present-day metal-rich globular clusters. 



In Sects. 12.11 and 12.21 we summarise our model, while the initial 
conditions of the simulations are presented in Sect. 12.31 The evolu- 
tion of the star cluster population in galaxy mergers is assessed in 
Sect. [3] where we also address their sensitivity to model parame- 
ters. We end this paper with a summary of our conclusions. 



2 SUMMARY OF THE MODEL 

We model the formation and evolution of star clusters semi- 
analytically, co upled to a numeri cal simulation code for galaxy evo- 
lution (STARS, Pelupessv 2005). Here we provid e a summary of 
the mod el, which was presented and validated by Kruiiss en et al.l 
l l2011bh . 
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2.1 Galaxy evolution and star cluster formation 

The evolution of the stellar and dark matter components are gov- 
erned by pure collisionless New tonian dynamics, ca lculated us- 
ing the Barnes-Hut tree method dBarnes &Hutlll986t) . The parti- 
cles sample the underlying phase space distribution of positions 
and velocities and are smoothed on length-scales of approximately 
0.2 kpc to maintain the collisionless dynamics and to reduce the 
noise in the tidal field (which is used for the cluster evolution, see 
Sect. 12.2b . The Euler equations for the gas dynamics are solved 
using smoothed particle hydrodynamics, a Galilean invariant Lan- 
grangian method f or hydrodynamics based on a particle represen- 
tation of the fluid (Monaghan 19921) . in the conservative formula- 
tion of ISpringel & Hemquistl ( 20021) . This is supplemented with a 
model for the thermodynamic evolution of the gas in order to rep- 
resent the physics of the interstellar medium (ISM). Photo-electric 
heating by UV radiation from young stars is included (assuming 
optically thin transport of non-ionizing photons). The UV field is 
calculated from stellar U V luminosities derived f rom stellar pop- 
ulation synthesis models ( iBruzual & Charlotll2003l) . Line cooling 
from eight elements (the main constituents of the ISM H and He 
as well as the elements C, N, O, Ne, Si and Fe) is included. We 
calculate ionization equilibrium including cosmic ray ionization. 
Furthe r de tails of the ISM m odel can be found in IPelupessv et al.l 
d2004l) and lPelupessvl d2005h . 

Star formation is implemented by using a gravitational insta- 
bility criterion based on the local Jeans mass Mj: 



2\V2 



< M re 



(1) 



where p is the local density, s the local sound speed, G the gravi- 
tational constant and M re t a reference mass-scale (chosen to corre- 
spond to observed giant molecular clouds). This selects cold, dense 
regions for star formation, which then form stars on a time-scale r s f 
that is set to scale with the local free fall time tg: 



Tsf = fsftff = 



(2) 



where the delay factor f s f « 10. Numerically, the code stochasti- 
cally spawns stellar particles from gas particles that are unstable 
according to Eq.Q]with a probability 1 — exp (— dt/r s f). The code 
also assigns a formation time for use by the stellar evolution li- 
brary, and sets the initial stellar and cluster population mass distri- 
butions (see below). Mechanical heating of the interstellar medium 
by stellar winds from young star s and supernovae is implemented 
by m eans of pressure particles ( Pelupessy et al. 2004]; IPelupessvl 
2005), which ensures the strength of feedback is insensitive to nu- 
merical resolution effects. In this way, the global efficiency of star 
formation is determined by the number of young stars needed to 
quench star formation by UV and supernova heating, which is set 
by the cooling properties of the gas and the energy input from the 
stars. 

For the purpose of this paper, in which the formation rate 
of star clusters is compared to their destruction rate, it is essen- 
tial that we obtain reliable estimates of the star formation rate 

(SFR). Our model for star formation reproduces the Kennicutt- 

l ll ~i i ii ~~i 

Schmidt (Schmidt 1959; Kennicutt 1989) pattern of star formation 

dGerritsen & Ickd 1 19971 ; IPelupessv 2005) and also gives realistic 
represen tation of the relation between molecular H2 and star for - 
mation dPelupessv et alj[2006l ; IPelupessv & Papadopoulosl 1200^) . 
On the other hand our model does simplify the star formation pro- 
cess considerably and this should be kept in mind. First, the ab- 



solute scaling of the star formation rates is somewhat uncertain 
and depends on the choice of parameters. It is mainly sensitive 
to effective feedback strength, but the feedback strength parame- 
ter has been independently constrained within a factor of two by 
considering the power spectra of the res ulting HI distribution maps 
dPelupessv et alj|2004l ; IPelup essv 2005). Secondly, the Jeans mass 
argument we use for our star formation model suffers from some 
limitations. Apart from the fact we take a single reference cloud 
mass M le f (in principle a more sophisticated model using a cloud 
spectrum could be constructed), it is also a strictly local criterion: 
this means that a given point in our simulation is either star forming 
or not regardless of the immediate environment. A more realistic 
star formation criterion would try to identify GMC-like structures 
in the gas distribution and then convert these into stars using its bulk 
properties (like mass, radius, irradiation and angular momentum) - 
possibly on the basis of more detailed modelling results of single 
GMC calculations. Lastly, our star formation model is based on 
the presence of a two phase interstellar medium - and while com- 
parison to structural properties of actual galaxies gives good re sults 
dPelupessv et alj2004l200r3 ; |Pelupessv & Papadopoulosl2009l) our 
use of smoothed-particle hydrodynamics (SPH) has known limita- 
tions in the representation of stro ng shocks and ins tabilities such as 
the Kevin-Helmholtz instability l Agertz et al .120071) . which are im- 
portant in two phase media. This could be checked in future work 

with shock resolving adaptive mesh refinement (AMR) or moving 

I ■ II 

mesh methods (such as AREPO, Springel 2010). 

Whenever a new star particle is spawned, a 'sub-grid' set of 

star clusters is generated. Their masses are drawn from a power 

law ICMF with an exponential truncation dSchechterlll976l) : 



NdM oc M exp(-M/M*)dM, 



(3) 



where N is the number of clusters, M is the cluster mass, and M* = 
2.5 x 10 6 Mq is the exponential truncation mass, which is con- 
sistent with high- mass end of the present day mass distribution of 
globular clusters dFail & Zha ng 2001 ; Kruiissen & Portegies ZwarJ 
120091) . T his reflects the observed mass distribution of young star 
clusters dzhang & Falll 1 19991 ; lLada & Ladal 120031 ; iLarsenl 120091 : 
Portegi es Zwart et al hoioh and is likely also the ICM F of the ma- 
jority of globular clusters dKruiissen & Coopeill201lh . We adopt a 
minimum cluster mass of Mmi n = 10 2 M . The cluster formation 
rate is assumed to be proportional to the SFR by adopting a con- 
stant cluster formation efficiency (CFE) of 90%, which is chosen 
to minimise Poisson noise. Because it is taken to be constant, the 
precise value of the CFE is irrelevant and acts as a normalisation 
of the number of clusters. The remaining 10% of the mass is con- 
sidered to be formed as unbound associations or field stars. In the 
simulation, the field stars are not physically separated from the star 
clusters, as each star particle contains both clusters and field stars. 
Because our cluster model is sub-grid, we presently cannot include 
clusters more massive than about 10 59 Mq, which corresponds to 
the adopted particle mass (see Sect. 12. 3t . The number of particles 
in the simulation was chosen to cover the cluster mass range of 
interest, while ensuring sufficient numerical resolution. 



2.2 Star cluster disruption 

The further evolution of the stellar clusters is computed with th e 
SPACE cluster models (Kruiissen & Lamers 2008; Kruiissen 2009), 
which include a semi-analytical description of the evolution of the 
cluster mass and its stellar conte nt. SPACE includes stellar evolu- 
tion from the Padova isochrones dMarigo et ali f2008). stellar rem- 
nant production, remnant kick velocities, dynamical disruption and 
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the evolution of the stellar mass function within the cluster due to 
the stellar mass dependence of the escape rate. The cluster evo- 
lution model has been coupled to properties of the tidal field by 
iKruiissen et al, (2011b) to include tidal evaporation and heating by 
tidal shocks. 

After their formation, the mass evolution of individual clusters 
is governed by mass loss due to stellar evolution and dynamical 
disruption: 

fdM\ /dM\ /dM\ 

with M the cluster mass and the subscripts 'se' and 'dis' denot- 
ing stellar evolution and disruption, respectively. The mass loss 
due to stellar evolution is obtained by taking the decrease of the 
maximum stellar m ass over one time step from the Padova models 
dMarigo et ai1l2008h . and integrating the mass function within the 
cluster over the corresponding mass interval. Upon the removal of 
these massive stars, the masses of their stellar remnants are added 
to the cluster mass. The dynamical mass loss is caused by two si- 
multaneous mechanisms. Firstly, the stars in the cl uster are drive n 
over the tidal boundary due to two-body relaxation dSpitzerll 19871) . 
Secondly, stars can gain energy from tidal shocks, i.e. fluctuations 
of the tidal field caused by passages through dense gal actic regions 
such as gian t molecular clouds (GMCs) or spiral arms (Gieles et al. 
I200dl2007l) . 

We parametrize the mass loss due to disruption as 

/dM\ _ /dM\ /dMA M M 

where 'rlx' and 'sh' denote two-body relaxation and tidal shocks, 
f d!s represents the time-scale for disruption by two-body relax- 
ation, and tfr S the time-scale for disruption by tidal shocks. Both 
time-scales are related to the tidal field. The derivation is given in 
IKruiissen et al.l d2011bl) . but here we give the final expressions. For 
t d j* the expression is: 

= 1.7 Gyr ( ^] , (6) 

4 V10 4 Gyr- 2 y 

where M4 is the cluster mass in units of 10 4 Mq, 7 = 0.62 is 
the ma ss dependence of the disruption time-scale dLamers et al.l 
l2005al) . which has a weak d ependence on the density profile of the 
cluster dLamers et al. 2010), and T is the tidal field strength. The 
tidal field strength is taken to be the largest eigenvalue of the tidal 
tensorQ which is determined by numerical differentiation of the 
force field. Gravity is smoothed on a length-scale of 0.2 kpc and the 
differentiation interval is 1 per cent of the smoothing length, which 
ensures that the influence of discret eness noise on cluster dis ruption 
is negligible. This was illustrated in lKruijssen et alJd2011bl) . where 
we also showed that the disruption of clusters due to tidal evapora- 
tion and tidal shocks is unaffected by passages of single particles 
for our choice of smoothing length and particle mass. Instead, their 

2 By doing so, we ignore potential second-order effects due to the other 
eigenva lues and the time evolution of the direction of th e largest eigen- 
vector I Tanikawa & Fukushige 2010; Renaud et al. 2011). This choice is 
made because the erratic tides in galaxy mergers with a gas component 
obstruct a straightforward implementation of these effects. None the less, 
the influence on our results should be minor for two reasons. Firstly, our 
model gives good agreement with di rect N-body simulations of cluster evo- 
lution iBaumgardt & Makino 2003). Secondly, the vast majority of clus- 
ter disruption is due t o tidal shocks instead of the steady tidal field (see 
IKruiissen et alj201 lbl and Sect.l3~2l. 



disruption is governed by the tidal influence of structures that are 
well-resolved with our resolution (also see Sect. l2.3t . The resulting 
small influence of numerical resolution on our results is verified in 
Sect. 1331 

For th e disruption t ime- sc ale due to tidal shocks, t he ap- 
proaches of IGieles et alj d2007l) and IPrieto & Gnedirj d2008l) can 
be combined to obtain dKruiissen et alj201 lbl) : 

«£ = 3.1 Gyr M 4 V7 — V ' ( to) , (7 ) 

where rh is the half-mass radius, J t id is the tidal heating parame- 
ter (seelGnedin et al.ll 19991 : IPrieto & G nedin 20081: IKruiissen et al.l 
l2011bl) . which follows from the integration of the tidal field over 
the duration of a shock, and At the time since the last shock. It 
reflects the time-scale on which the cluster is heated and is deter- 
mined individually for each component of the tidal tensor by identi- 
fying local minima with s ufficient (lg) contrast with respect to the 
preceding maximum (see lKruiissen et al]|2011bl) . Because the dis- 
ruption time-scale due to tidal shocks depends on cluster density, 
it is important to includ e a descripti o n for the half-mass radius. It 
was recently shown by IGieles et al.l feoill) that cluster radii pass 
through two evolutionary phases. Initially, a cluster expands to fill 
its tidal boundary, during which time the half-mass relaxation time 
remains constant, i.e. rh oc M -1 ^ 3 . After filling its tidal boundary, 
the cluster continues in the 'mass-loss dominated regime' along 
tracks of rh oc M x , with x = 1/6 to 1/3 depending on the es- 
cape criterion. The duration of the first phase depends on the initial 
conditions of cluster formation, while the second phase lasts un- 
til the total disruption of the cluster. Since our models have been 
tu ned to agree with the JV-body simulations of cluster disruption 
by Baumgardt & Makino d2003l) , we assume an evolution of the 
half-mass radius r h = 4. 35 pc (M/10 4 M fl ) om , which is consis- 
tent with their work (see lKruiissen et al]|201 lbl). This rel ation lies 
in the second evo lutionary phase from | Gieles et alj d201ll) . because 
the clusters from IBaumgardt & M akino (2003) are initially filling 
their tidal boundaries. Because the initial conditions of cluster for- 
mation and their impact on the mass-radius relation are quite un- 
certain, we validate our results using other mass-radius relations in 
Sect. |331 

Both Eqs. [6] and [7] are calibrated for clusters with a King pa- 
rameter of Wo = 5. For other density profiles, the constants in the 
equations change, but the lifetimes of the clusters are similar. They 
have been compared and c alibrated to the JV-bo d y sim ulations of 
star cluster disr uption by [B aumgardt & Makino (2003) to ensure 
their accuracy dKruiissen et al.l 1201 lbl) . While most of their sim- 
ulations concern clusters on circular orbits, we have used those 
simulations of clusters on eccentric orbits to verify our models for 
changing tidal fields. 

Like the cluster formation rate, it is essential for the purpose 
of this paper that the estimated disruption rates are reliable. In 
IKruiissen et ai] d2011bh . we have therefore tested our model for a 
range of different cosmic settings including several isolated disc 
galaxies and different kinds of galaxy mergers, and compared the 
results to observations. The simulated age distributions of star clus- 
ters in disc galaxies and their correlation with gala ctocentric ra- 
dius ar e in accordance with observational results (see Bastian et al.1 
2011a for an analysis of the cluster population of M83). We also 
found that in isolated disc galaxies with 15-30% of their baryonic 
mass in gas, typically 85% of the cluster disruption is accounted 
for by tidal shocks (Eq.[7), while the remainder is covered by two- 
bo dy relaxation (Eq.lp. Th is is in excellent agreement with a study 
bv lLamers & Gielesl (l2006), who found that in the solar neighbour- 
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Table 1. Details of the initial conditions for the disc galaxy models. 



ID 






z 


A 






jystar 

□ISC 


N star 

bulge 


M halo a 
-'"part 


M baryi ' 


Comments 


IdA 


0.20 


10 12 


2 


0.05 


10 6 


10250 


41000 


10000 


10 6 


8 x 10 5 


low gas fraction 


ldB 


0.30 


10 12 


2 


0.05 


10 b 


15375 


35875 


10000 


10'' 


8 x 10 5 


standard model 


ldC 


0.50 


10 12 


2 


0.05 


10 6 


25625 


25625 


10000 


10 6 


8 x 10 5 


high gas fraction 


ldD 


0.30 


5 x 10" 


2 


0.05 


5 x 10 s 


7688 


17938 


5000 


10 6 


8 x 10 5 


half mass 


IdE 


0.30 


10 12 


2 


0.05 


10 6 


15375 


35875 





10 6 


8 x 10 s 


no bulge 


ldF 


0.30 


10 11 


2 


0.05 


10 6 


15375 


35875 


10000 


10 s 


8 x 10 4 


low mass 


ldG 


0.30 


10 12 


2 


0.10 


10 6 


15375 


35875 


10000 


10'' 


8 x 10 5 


high spin 


ldH 


0.30 


10 12 





0.05 


10 6 


15375 


35875 


10000 


10 6 


8 x 10 5 


low concentration 


ldl 


0.30 


10 12 


5 


0.05 


10 6 


15375 


35875 


10000 


10 6 


8 x 10 5 


high concentration 



a In solar masses (Mq). 



hood about 80% of the disruption is contributed by tidal shocks. 
The high relative contribution of tidal shocks to cluster disruption 
shows that the tidal field in an isolated galaxy is far from smooth 
due to encounters with GMCs and spiral arms. This is an important 
similarity to galaxy mergers, in which the tidal field also varies, 
albeit to a larger extent. We have also applied our models to the 
Antennae galaxies (Kruijssen & Bastian, in prep.) and find good 
agree ment with the observ ed cluster age and mass distributions 
from lWhitmore et"al] (2007). These results provide a good starting 
point to apply our model to galaxy mergers and follow the forma- 
tion and evolution of the entire star cluster population for different 
galactic histories. 



2.3 Initial conditions 

We use the set of simulations described in Kruijssen e t alj d201 lbl) 
and summarise the adopted parameter sets here. The simulations 
follow the evolution of the star cluster population in isolated disc 
galaxies and galaxy mergers. The disc galaxies are generated with 
parameters that can be related to the outcomes of cosmological 
ACDM galaxy formation models dMo et al.lll998l: ISpringel et alj 
120051) , They consist of a dark halo with a lHernquistl dl990h profile, 
an exponential stellar disc, a stellar bulge (except for one model) 
and a thin gaseous disc, construct ed to be in self gravit ating equi- 
librium if evolved autonomously dSpringel et ail 12005). The dark 
matter haloes have concentration parameters related to their to- 
tal ma sses and condensation redshifts according to Bullo ck et alj 
d200lh . implying that for a fixed mass the halo concentration in- 
creases with redshift. The total mass is related to the virial ve- 
locity Vvir and the Hubble constant H(z) at redshift z as Mm = 
V^ n /[l0GH(z)]. For all galaxies, the baryonic disc is constituted 
by a gaseous and stellar component, having a mass fraction im = 
0.041 of the total mass, while the bulge (when included) consists 
of a stellar component only, having a mass fraction nib = 0.008 
of the total mass. The fraction of total angular momentum that 
is constituted by the disc (jd) is taken identical to m&. The scale- 
length of the bulge and the vertical scale-length of the disc are 0.2 
times the radial scale-l ength of the dis c, which is determined by 
the degree of rotation dMo et al.|[l998h through the spin parame- 



Table 2. Details of the initial conditions for the galaxy merger models. 



ter A = J\E\/GM^ 2 , in which J is the angular momentum of the 
halo and E its total energy. We have chosen the parameter sets to 
cover a reasonable spread in galaxy properties, specifically the gas 
fraction of the baryonic disc f gil!i , their total mass M v h , the spin pa- 
rameter A and the presence of a bulge. The resulting disc galaxy 
model parameters can be found in TableQ] which lists f gas , M v ; r , A, 
the number of particles in the different components of the model 



ID 


Discs 


01 


</>! 


8 2 


02 




Type b 


1ml 


AA 
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PP 


lm2 


BB 














6 


PP 


lm3 


CC 














6 


PP 


lm4 


BD 














6 


PP 


lm5 


EE 














6 


PP 


lm6 


FF 














6 


PP 


lm7 


GG 














6 


PP 


lm8 


HH 














6 


PP 


lm9 


II 














6 


PP 


lmlO 


BB 


60 


45 


-45 


-30 


6 


PPi 


lmll 


BB 


180 











6 


PR 


lml2 


BB 


-120 


45 


-45 


-30 


6 


PR, 


lml3 


BB 


180 





180 





6 


RR 


lml4 


BB 


-120 


45 


135 


-30 


6 


RRj 


lml5 


BB 














12 


PPw 


lml6 


CG 


-120 


45 


-45 


-30 


10 


PRi 


lml7 


BB 








71 


30 


6 


Barnes 


lml8 


BB 


-109 


90 


71 


90 


6 


Barnes 


lml9 


BB 


-109 


-30 


71 


-30 


6 


Barnes 


lm20 


BB 


-109 


30 


180 





6 


Barnes 


lm21 


BB 








71 


90 


6 


Barnes 


lm22 


BB 


-109 


-30 


71 


30 


6 


Barnes 


lm23 


BB 


-109 


30 


71 


-30 


6 


Barnes 


lm24 


BB 


-109 


90 


180 





6 


Barnes 



a In kpc. All angles are in degrees. 

''Indicates prograde-prograde (PP), prograde-retrograde (PR), 
retrograde-retrograde (RR) or 'Barnes' (see text). Subscripts 'i' and 
'w' denote inclined/near-polar and wide orbits, respectively. 



galaxies, and the particle masses of the halo particles Mpa" 1 , and 
baryonic particles . For each set of initial conditions, M^" 
and Mpa, J t y are chosen to be very similar. Our particle resolution is a 
trade-off between enabling the formation of high-mass star clusters 
(which are limited by the particle mass due to their sub-grid treat- 
m ent, see Se ct.[2T} and resolving the galaxy dynamics. We verified 
in lKruiissen et alj d2011bl) that the adopted resolution is sufficient 
to reliably model the cluster disruption, because the influence of 
encounters with single particles is negligible compared to the tidal 
perturbation of clusters by more massive, resolved structures. 

In the galaxy merger simulations, the disc galaxies follow Ke- 
plerian parabolic orbital trajectories with initial separations of ap- 
proximately 200 kpc. The actual orbit will decay due to dynamical 
friction, which leads to the merging of the galaxies. The orbital ge- 
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Figure 1. Evolution of the star cluster population during a galaxy merger (simulation 1ml 1 from Table [2). The surface density of the gas is displayed in 
greyscale, while the particles that contain star clusters are shown in colours denoting the ages of the clusters as indicated by the legend. The subsequent 
images show the collision at eight characteristic moments: t = 0. 1 Gyr, briefly before the first passage; t = 0.3 Gyr, just after the first passage; t = 0.8 Gyr, 
in between the first and second passage; t = 1.5 Gyr, just before the second passage; t = 1.6 Gyr, just after the second passage; t = 1.7 Gyr just before 
the final coalescence; t = 2.1 Gyr, during the coalescence and the infall of remaining gas clouds; t = 3.8 Gyr, when only a merger remnant is left. See 
http : //www. mpa-garching.mpg.de/-diederik/lmll clusters .html for a movie of the full time sequence. 



ometry of an interaction is characterised by the directions of the 
angular momentum vectors of the two galaxy discs and the peri- 
centre distance of the parabolic orbit J? pcl i. The angular momentum 
vectors of the galaxies are determined in spherical coordinates by 
angles 6 (rotation perpendicular to the orbital plane) and (f> (rota- 
tion in the orbital plane). These and other relevant parameters are 
listed in Table [2] which covers three subsets of simulations. The 
first set of eight runs follow a common configuration, in which the 
discs rotate in the orbital plane. They are used to test the influence 
of the properties of the progenitor discs. The six subsequent runs 
are aimed at investigating the impact of orbital parameters on the 
star cluster population. We rotate the progenitor discs to include 
retrograde rotation and near-polar orbits, which represent the op- 
posite extreme with respect to the co-planar configurations of the 
first eight runs. A wider orbit and a 'random' major merger are also 
considered. The third group con tains the eight 'ran dom' configura- 
tions from lHopkins et ail J2009h (see lBarnesl 19881) . which together 
sample the phase space of possible orbital geometries. 

All galaxies are generated without any star clusters, and we set 
t = after 300 Myr of evolution to initialise the cluster population. 
As described in Sect. El the cluste r s have masses between 10 2 and ~ 
10 59 M Q , following a lSchechterl ( fl97rj) type initial mass function. 
The chemical composition of the clusters is set to solar metallicityd 
and we assume a King parameter of Wo = 5. 



3 This choice only affects the stellar evolutionary mass loss. 



3 EVOLUTION OF THE STAR CLUSTER POPULATION 

In this section, we apply our model to the evolution of the star clus- 
ter population in galaxy mergers. We start out by explaining an 
illustrative example, before discussing the other simulations and 
trends with merger properties. 



3.1 Illustrative example 

Figure [TJ shows a classical sequence of the evolution of a galaxy 
merger simulation together with the results from our cluster evo- 
lution model (simulation 1ml 1 from Table|2j. The panels in Fig.[TJ 
show the distributions of gas and star clusters at different times dur- 
ing the interaction. The first image displays the galaxies as they ap- 
proach each other for their first passage (at t = 0.1 Gyr), when the 
tidal interaction between the galaxies is still relatively weak and the 
SFR is at a low-to-intermediate level (~ 6 Mq yr~'). The spatial 
distribution of star clusters is restricted to both galaxy discs, where 
the gas resides from which they are formed, and their destruction 
rate is still low since it is only driven by the internal galactic tidal 
field and encounters between clusters and GMCs. 

In the second and third images of Fig. [T](t = 0.3-0.8 Gyr), 
the galaxies are shown (briefly) after their first passage. By this 
time, the gravitational interaction has produced extended tidal tails. 
Most star clusters still follow the morphology of the gas because 
they have just been formed in a large starburst (about 50 Mq yr _ ) 
that was triggered by the angular momentum loss and consequent 
inflow of the gas during the first pericentre passage. In the second 
image, the total number of clusters reaches a peak, with an increase 
of ~ 40% with respect to the first panel, but in the third image 
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only one third of this number is This is caused by the large 
central gas density that drives the starburst, prompting a stronger 
increase of the tidal perturbation of star clusters than of their for- 
mation rate. Some intermediate age clusters have been ejected from 
the discs by the interaction. They represent the first star cluster con- 
stituents of a stellar halo forming around the two galaxies. The 
mechanisms of cluster migration an d natural selection that were 
identified in lKruiissen et alj ( feOllbl) are evident here: clusters are 
escaping the dense regions in which they were formed and those 
clusters in quiescent, low-density parts of the galaxies have higher 
chances of survival. As a result, the mean disruption rate decreases 
with cluster age. The enhanced disruption of young clusters due to 
their tidal interaction wit h the primordial environm ent was named 
the cruel cradle effect in iKruiissen et al.1 d201 lal) . and can affect 
clusters up to ages of r ~ 200 Myr, contrary to the early disruption 
of clusters by gas expulsion ('infant mortality' ), which takes place 
on a gas expulsion time-scale of ~ 10 Myr dLada & Ladall2003l ; 
iGoodwin & Bastiarl2006l|Pelupessv & Portegies Zwartll201lh ! 

As the galaxies proceed to merge, the effects of the interaction 
intensify. The fifth and sixth panels of Fig. Q] display the galaxies 
during the short interval between the second passage and their final 
coalescence (t = 1.6-1.7 Gyr), in a confi guration that is s imilar to 
the Antennae' galaxies (NGC 4038/9, see lKarl et alj20ld) . During 
this phase, the remaining gas is funnelled towards the centres of the 
galaxies, where it cools to form large numbers of stars and star clus- 
ters. This second starburst is accompanied by an even stronger in- 
crease of the cluster destruction rate, this time disrupting well over 
50% of all clusters. A large number of clusters is ejected from the 
central region into the stellar halo that surrounds the galaxies. Away 
from the turmoil, these clusters will be able to survive the coales- 
cence of the galaxies. By the time the merger is completed, most of 
the surviving clusters will have form ed at this moment or a round 
the time of the first snapshot in Fig.fTl fKruiissen et al.ll201 fbT) . 

When the merger is completed, as is shown in the last image of 
Fig.[TJ(t = 3.8 Gyr), the system has transformed into a giant ellip- 
tical galaxy, in which the star cluster system has dispersed into the 
stellar halo. The formation rate of clusters drops to a minimum af- 
ter the merger, due to the depletion of the gas during the starbursts. 
The gas depletion also affects the disruption rate: it causes the typ- 
ical tidal shock strength to decrease after the coalescence of the 
two galaxies, implying that massive, dense clusters are no longer 
affected by disruption. However, the migration of clusters towards 
radial orbits during violent relaxation causes the tidal shock heat- 
ing to remain high, as clusters on very eccentric orbits are being 
disrupted by bulge shocks. At this stage, low-mass star clusters are 
preferentially disrupted, and their destruction leads to an increase 
of the mean cluster mass. 



3.2 Formation versus destruction 

We have tested the generality of these results by analysing the full 
set of simulations from Sect. 12.31 The results of two simulations are 
shown in Fig. [2] where the star formation history (SFH) as well as 
the time-evolution of the mean tidal shock heating, the mean am- 
bient densities of gas and starsQ and the number of star clusters 

4 The total number of clusters depends on the lower mass limit (see 
Sects. I3~2l and l3.3t . which in our simulations is taken to b e 10 2 Mpy 

5 The se are determined by using the approach from ICasertano & Hull 
1 1985), where the density is averaged over a sphere with radius equal to 
the distance to the Nth nearest neighbour. We adopt N = 7. 
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Figure 2. Evolution of the star formation rate (SFR, first row), mean tidal 
heating ((itid). second row), mean ambient gas, stellar and total densi- 
ties ((pamb), third row), and the number of star clusters (N, fourth row) 
for two different galaxy merger simulations. The left-hand panels show 
the results for the prograde-retrograde encounter from Fig. [T] (simulation 
lmll from Table [2), while the right-hand panels represent a similar en- 
counter with both galaxies rotating in the prograde direction (anticlock- 
wise in the configuration of Fig. [T] simulation lm2 from Table |5J. The 
thick dots mark the moments that are displayed in Fig. [T] The number of 
star clusters in the bottom panels is shown for different cluster mass cuts 
(log(M/M0) > {2,3,4,5}). The bottom lines also include an age limit 
(t < 10 Myr), and show the number of young massive clusters as a func- 
tion of time. The dotted curves denote the results for the two disc galaxies 
evolving in isolation. The vertical dashed lines indicate the times of first 
and second pericentre passage and the shaded areas specify the time inter- 
val over which the final coalescence occurs. 

for different cluster mass cuts are shown. The figure also includes 
a comparison with the two disc galaxies evolving in isolation. Just 
after the pericentre passages, the galaxies exhibit a pronounced in- 
crease of the SFR (0.5-1 dex), but an even stronger increase of the 
mean tidal shock heating (1—1.5 dex), leading to a decrease of the 
total number of clusters by nearly two orders of magnitude towards 
the end of the simulations. It is not the tidal field strength itself 

6 In fact, the mean tidal field strength is lower in the merger simulations 
than in the isolated disc galaxies. This is caused by several factors, but is 
mainly related to the rapid destruction of clusters by tidal shocks, which 
works most efficiently in the regions of a merger where the absolute tidal 
field strength is also higher. The resulting cluster population is biased to- 
wards larger galactocentric radii than in the isolated disc case, implying 
that the mean tidal field strength is lower. Being the dominant source of 
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but the frequency and strength of tidal shocks which leads to the 
enhanced disruption of clusters. 

If the (environmentally dependent) increase of the tidal dis- 
ruption is neglected, the total number of clusters increases by a fac- 
tor of 5-6 during the first pericentre passage for simulation 1ml 1 
compared to the progenitor discs evolving in isolation. For both 
galaxy mergers in Fig. [2] the destruction of clusters is most promi- 
nent after the second pericentre passage and during the coalescence 
of the galaxies. As indicated in Sect. 13. II the decrease of the num- 
ber of clusters is largest for the lowest cluster masses, which is 
clearly seen in the number evolution for different cluster mass cuts 
in Fig. [2] After the mergers are completed, their SFRs become 
lower than would have been the case had the galaxies evolved in 
isolation. At this stage, the mean tidal shock heating remains rel- 
atively high due to the disruption of clusters on radial orbits (see 
below). 

The densities in Fig. [2] illustrate that both the starbursts and 
the episodic increase of the tidal shock heating are caused by the 
growth of the ambient gas density. The stellar density is almost al- 
ways higher than the gas density, but during the pericentre passages 
and final coalescence, the development of peaks in (Jtia) correlates 
with (pgas) rather than (p s tar): contrary to (p s tar), which remains at 
a constant, high level after each increase, (p ga s) and (Im) return to 
lower values. The reason that the tidal shock heating is dominated 
by the ambient gas density instead of the (higher) stellar density is 
that the gas is more structured, which produces faster tidal shocks 
that cannot be absorbed by the adiabatic expansion of the clusters. 
Only after the final coalescence of the galaxies, the tidal shock heat- 
ing is dominated by the stellar density. This occurs after t = 3 Gyr, 
when most of the cluster disruption is caused by the tidal shocking 
of clusters on highly eccentric orbits, which are falling in from the 
tidal tails or have migrated to a radial orbit due to violent relaxation. 

In time, the secular evolution of a merger remnant decreases 
the orbital anisotropy of the surviving clusters through the dis- 
ruption of clusters on eccentric orbits. To quantify the orbital 
anisotropy of the star cluster system as a function of galactocen- 
tric radius i? gc , the anisotropy parameter is defined as 



/3(R g c) = 1 



2<v, 2 )' 



(8) 



where (v r 2 ) is the mean square radial velocity in a radial bin cen- 
tered at J? gc , and (v t 2 ) is the mean square tangential velocity. For 
the isotropic case we have f3 = 0, while /3 > and j3 < indicate 
preferentially radial and tangential orbits, respectively. In Fig. [3] 
the time evolution of the anisotropy parameter is shown as a func- 
tion of galactocentric radius, for snapshots during and after the final 
coalescence in simulation lml 1. Up to t = 3.8 Gyr, the kinematics 
of the inner ~ 10 kpc are dominated by rotation because the en- 
counter is co-planar. Outside this radius, the cluster system quickly 
becomes radially anisotropic, with an anisotropy radius close to 
R A = 30 kpc. As mentioned earlier, this is caused by the infall 
of clusters from the tidal tails and the migration to eccentric or- 
bits due to violent relaxation. However, after t = 3.8 Gyr the radial 
anisotropy disappears due to the destruction of clusters on radial 
orbits (also see Fig. |2j- In the inner ~ 10 kpc, the cluster system 
also evolves towards isotropy. The anisotropy radius increases to 
Ra = 50-60 kpc, similar to t he result found fo r globular clusters 
from accreted dwarf galaxies jPrieto & Gne din 2008). This shows 



the cluster destruction, the mean tidal shock heating is not affected by this 
selection effect (see Fig.ff). 




t = 2.5 Gyr 
t = 3.8 Gyr 
t = 4.9 Gyr 



10 1 
R g c [kpc] 

Figure 3. Orbital anisotropy of the surviving star clusters ft (see Eq. [8) 
as a function of galactocentric radius. The relation is shown at different 
times (indicated by the legend) during the final coalescence and in the 
merger remnant of simulation lmll. The clusters are binned using an equal 
numbers of clusters per bin, with the error bars denoting the error on the 
mean. Preferentially radial and tangential orbits are separated by the hor- 
izontal dashed line at = 0, which indicates orbital isotropy. The dotted 
line is included for reference and shows the radial dependen ce of /3 for 
the parametrization of the velocity ellipsoid from Aguilar et al. 1 1988) with 
anisotropy radius Ra = 30 kpc. There are indicatio ns that the globular clus- 
ter system of M87 has a similar anisotropy profile iStrader et alj201 H) . 



that it is not straightforward to distinguish between in-situ and ex- 
situ formation based on the orbital (an)isotropy of a cluster system. 

The results in Fig. [2] suggest that galaxy mergers efficiently 
disrupt star clusters, in apparent contradiction with the rich star 
cluster population s that are observed in colliding galaxies such 
as th e Antennae dWhitmore et al] fl99^) and M51 teastian et al.l 
l2005h . However, our above analysis concerns the entire star clus- 
ter population in a merger, while observations are naturally con- 
strained to bright clusters, which are typically massive and young. 
When limiting our results to the young (< 10 Myr) and massive 
(> 1 4 -10 5 Mq) cluster s that are easily detected in observations 
(e.g. IZhang & Falll 1 19990 . Fig. [2] shows that the number of clus- 
ters that would be 'observed' from our simulation temporarily in- 
creases by more than a factor of three during starbursts. The in- 
crease i n observed galaxy mergers may be even higher than thi s 
(see e.g. lSchweizer et alJl996r , lMiller et alJl997r . lZepf et al.ll999T) . 
but the factor of three increase is a lower limit for a number of rea- 
sons. Firstly, simulation lmll is one of the more monotonously 
disruptive merger simulations in our sample (second from the right 
in Fig. [4] see below). It can be contrasted with simulation lm7, 
which exhibits the highest degree of variation over the course of 
the merger: the number of young massive clusters is temporar- 
ily boosted by a factor of 6-15 during the time interval t = 0.5- 
2.2 Gyr, whereas the total number of clusters eventually settles at 
only 7% of the amount the progenitor galaxies would have had in 
isolation. Secondly, we did not include a vari able CFE, which may 
increase with the s tar formation rate density dGoddard et al 1 l20ld 
lAdamoetai]|201ll. although see Isilva- Villa & Larserj|201 lh . This 
could imply that the number of clusters increases by an additional 
factor of 2-3 during starbursts. Thirdly, our cluster masses are lim- 
ited to 10 5 ' 9 Mq, which obstructs the formation of the extremely 
massive (^ 10 6 Mq) clusters that are observed in galaxy mergers 
and for which the relative increase with respect to quiescent galax- 
ies is most evident. These effects could conspire to yield a transient 
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Figure 4. Star cluster 'survival fraction' f surv as a function of the starburst 
intensity parameter © (see text and Eq.|9). The corresponding logarithmic 
mean of the gas depletion time-scale (tiepl) = ©~'/ 2 of the two starbursts 
during a merger is indicated along the top axis. Symbols denote the 24 
merger simulations from Table [2] and the dotted line gives a power law 
fit. 



relative increase of young massive clusters during the starbursts of 
up to a factor of ~ 30. None the less, in terms of numbers, a star 
cluster population is dominated by the unseen low-mass star clus- 
ters that are effectively destroyed during a merger before they reach 
ages much older than a few tens of Myr. 

3.3 A generalised relation for star cluster survival 

For all simulations, the results are in accordance with those shown 
in Fig. [2] as they exhibit a very similar increase of the mean tidal 
shock heating and corresponding decrease of the number of clus- 
ters during the merger. The number of clusters in our merger rem- 
nants is always 2-50% of the amount that the two discs would have 
contained in isolation. Much of the variation is caused by the dif- 
ferent orbital geometries. Retrograde, co-planar encounters lead to 
enhanced angular momentum loss of the gas and correspondingly 
stronger starbursts and greater destruction of clusters, decreasing 
their number by up to a factor of 50. Galaxies on wide or inclined 
orbits such that they follow near-polar trajectories prompt a weaker 
effect due to a less pronounced gas inflow, yielding a decrease of 
about a factor of 2-10. 

We find that the total number of surviving clusters strongly de- 
creases with increasing peak SFR. This trend is a consequence of 
the disruptive power of dense, star-forming environments. To quan- 
tify this trend across all simulations, one can define the ratio of the 
number of clusters in the galaxy mergers relative to the number of 
clusters in the isolated progenitor galaxies (f smv ). We have tested 
the dependence of f slllv on several generalised forms of the peak 
SFR, by normalising SFR pea k to the galaxy stellar, gas or baryonic 
mass (Mgtar, M gas or Mbary, respectively). It is found that f surv most 
tightly correlates with SFRpeak /Mgas. The simulated galaxy merg- 
ers typically experience two starbursts (during the pericentre pas- 
sages), and therefore the value of f slllv in a merger remnant includes 
the effect of two starbursts. For the case of a single starburst, we 
write a power law formulation f smy = C[SFR P cak/(M gas yr -1 )]". 
For the number of clusters in a merger remnant this implies: 



fsurv — fsurv, 1 fsurv,2 — C 



SFR 



peak, 1 



SFRnea 



peak,2 



M gasa yr-' M gas-2 yr- 
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where subscripts 1 and 2 indicate the first and second starbursts, 
respectively, and we have defined a starburst intensity parameter 

EE SFRpeak.lSFRpeak^/CMgas.lMga^ yr" 2 ). 

Figure [4] shows f smv (0) when measuring f surv at t = 4.8 Gyr, 
which is typically 2 Gyr after the completion of each merger. The 
'survival fraction' f surv very clearly decreases with increasing star- 
burst intensity 0. A simple power law fit to the data points in Fig. [4] 
gives C = 4.5 ± 1.5 x 10 s and a = 0.79 ± 0.13 over almost two 
orders of magnitude in 0. While the deviation of some data points 
is as high as 0.5 dex, the correlation is quite remarkable consid- 
ering the wide range of boundary conditions that is covered (see 
Tables [TJ and [2](. The relation flattens when increasing the lower 
mass limit of the clusters, since massive clusters are less rapidly 
disrupted than low-mass clusters. For our models, this can be ap- 
proximated to reasonable accuracy by C = 4.5 x lO^M^j-j and 
a = -0.77 + 0.22 log M ra in.2, for M ral „,2 = M mm /10 2 M and 
10 2 ^ M m i„/M0 ^ 10 4 . At larger minimum masses, C and a 
remain constant, although it is uncertain to what extent this may 
be the result of our maximum mass limit of 10 59 Mq. Another 
source of uncertainty is the variation of the ICM F truncation mass 
M» with the galactic environme nt or SFR dBastianll200"& iLarsenl 
120091 iKruiissen & Coopeill201 ll) . If M* increases with the SFR, 
galaxy mergers naturally yield a net production of star clusters with 
masses larger than the value of M* in quiescent progenitor galax- 
ies. T his has been reported to be about M* ~2x 10 5 Mq ( ILarsenl 
2009), which thus indicates the mass scale that separates net clus- 
ter destruction at low cluster masses from a net production at higher 
masses. 

If we write Eq. [9] in terms of the gas depletion time-scale 
tdepi ee M gils /SFRpeak, the above results in a generalised expression 
for f sulv after a single starburst, which is given by 



f surv (M > M^) = 4.5 x 10~ 8 M- 



tdepl 

yr 



0.77-0.22 lo g M min 2 



,(10) 



c 2 e a , (9) 



for 10 2 < M mm /M Q < 10 4 and 0.1 < t dcp i/Gyr < 3. A naive 
extrapolation of this expression gives a net increase of the number 
of star clusters during starbursts above ~ 3 x 10 5 Mq, very sim- 
ilar to the approximate value of M* in quiescent galaxies. Since 
we neglect any environmental variation of M*, this similarity is 
a coincidence that potentially allows Eq. QJJ] to be extended to 
Mmin > 10 4 Mq. This will need to be addressed in a future work 
that does account for a variation of the truncation mass. Because 
t d ~ i is a measure of the intensity of the starburst, Eq. \W\ reflects 
that clusters are disrupted by the dense star-forming environment. 
Figure|4]has thus shown that, ironically, star formation kills. 

The enhanced cluster disruption rate during starbursts is also 
demonstrated by the displacement of the peaks in the cluster age 
distribution and SFH that was presented in lKruiissen et al. 
Figs. 15 and 17). Star clusters that are formed during starbursts are 
found to experience such an elevated disruption rate that they have 
severely lower survival chances than clusters formed in quiescent 
environments. As a result, the peaks in the cluster age distribution 
and SFH can differ by up to 200 Myr, with the bulk of the sur- 
viving star cl usters being formed pri or to the height of star forma- 
tion (also see lchien~& Barnes 2010). This implies that the cluster 
age distribution may indicate the occurrence of a starburst, but (de- 
pending on the strength of the starburst) cannot always be used to 
accurately determine its time or duration. For quiescent galaxies, 
the age distribution of star clusters does reflec t the SFH quite well , 
modulo a correction for cluster disruption (cf. lLamers et al.|[2005al : 
lBastianetalj2011bl) . 
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Figure 5. Evolution of the mass distribution of star clusters during merger 
simulation lml 1. Shown are the distributions at different times t. We refer 
to Fig.[T]for the merger episodes to which these times correspond. As time 
progresses, the distribution shifts downwards due to the net destruction of 
clusters. The slope of the initial mass distribution is shown as a dashed line, 
which would have closely resembled the mass distribution at all times had 
the two galaxies evolved in isolation. 

3.4 The cluster mass function 

The preferential destruction of the low-mass clusters causes the ini- 
tially scale-free (except for the Schechter-type truncation) cluster 
mass distribution to develop a characteristic massQ which is shown 
in Fig. [5]for simulation lmll. Most of this transformation occurs 
during the final coalescence of the galaxies from t = 1.8 Gyr on- 
wards, when the disruption rate is no longer high enough to affect 
the most massive clusters, but is still sufficient to efficiently destroy 
low-mass clusters. At t = 4.9 Gyr, the peak mass is about If/ 7 M 
and increases steadily. Tentative evidence for this peaked form of 
the cluster mass distributi on is also found in observatio ns of re- 
cent merger remnants (e.g. lGoudfrooii et al . 2004, 2007), albeit at 
higher masses (we refer the reader to Sect. |3.5| for a discussion of 
the variation of the modelled peak mass with model parameters). 

The surviving population of clusters that were formed before 
and during a merger bears hints of observed globular cluster sys- 
tems. First of all, the spatial configuration of these clusters is com- 
parable to th at of the globular cluster po pulation of the Milky Way 
( Harris! fl996h . gi ant elliptical galaxies fHarrisll200"9l) and young 
merger remnants dSchweizer et alj["l996l) . following a power law 
density profile with index —3.2 in the outer parts, which is the ap- 
proxi mate behaviour of a de Vaucouleurs profile Jde Vaucouleursl 
1948). Secondly, the development of a peak in the mass distribu- 
tion at a mass of 10 2 ' 7 Mq is very suggestive. This peak mass is 
still lower than the cha racteristic mass of globular cluster systems 
(10 5 MM. lHarris|[l99r3) . which would be attained during the sev- 
eral billions of years of star cluster disruption following a high- 
redsh ift merger until the present day JVesperinil200ll ; iFall & Zhand 
l200lMKruiissen & Portegies ZwarJ2009l) , possibly also due to sub- 
sequent collisions with other galaxies. Lastly, due to the high peak 

7 This requires that tidal shocks most efficiently disrupt low-mass clus- 
ters, i.e. that the density of clusters increases with their mass, and therefore 
would not occur for mass-radius relations rj, <x M* with <5 J: 1/3. How- 
ever, sucliasteongcorrekition is not supported b y observational evidence 
(e.g. Harris 1996; Larsen 2004; Bastian et al. 2005). We explore the depen- 
dence of the results on the mass-radius relation in Sect. 13. 51 
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Figure 6. Radial variation of the peak mass (top panel) and dispersion (bot- 
tom panel) of the cluster mass distribution at two different times in a merger 
remnant. Solid lines show the data from simulation lmll in bins with an 
equal number of clusters per bin, while dotted lines represent power law fits 
with slopes as indicated by the labels. 



SFR, a merger can produce a population of clusters that extends 
to higher masses than for isolated galaxies, in agree ment with ob- 
servations (Bastian 2008; Kruijssen & Coopej20l"lh . It is therefore 
capable of producing clusters with the initial masses needed to sur- 
vive for a Hubble time. 

Shortly after the completion of a merger, secular cluster dis- 
ruption increases the characteristic mass by 0.3-0.4 dex per Gyr for 
the next two gigayears (also see Fig.[6). If the merger took place in 
the early universe (> 9 Gyr ago), the characteristic mass would thus 
have the time to evolve to that of observed globular cluster systems. 
The combination of several globular cluster-like characteristics (the 
spatial distribution, characteristic mass, and high maximum cluster 
mass) would not be reproduced without the starburst and gas de- 
pletion, the migration of clusters into the halo, and the enhanced 
disruption occurring during the starburst. 

It is tempting to interpret the existence of a peak in the mass 
distribution of the surviving star clusters as the early formation of 
a globular cluster syste m. Such a scenario was first proposed by 
lAshman & Z epf ( 1992). If this were the case, the orbital kinemat- 
ics of the clusters should evolve towards a state in which there 
is no radial trend of the characteristic mass, as is the ca se for 
the gl obular cluster system s of the Milky Way (Harris 199^) and 
M87 dVesperini et alj |2003). This could proceed by orbital migra- 
tion or by the destruction of clusters with certain orbital charac- 
teristics. The violent relaxation occurring during galaxy mergers 
is ind eed efficient at ejecting clusters from their original environ- 
ment jPrieto & Gnedinl 120081 : iBastian et alj|2009l ; iKruijssen et al.l 
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Figure 7. Peak mass M pca k and logarithmic dispersion aw jvj as a function 
of the starburst intensity parameter © (left-hand panels) and of each other 
(right-hand panel). Symbols denote the 24 merger simulations from Table|2] 
The cluster mass functions are taken at t = 4.8 Gyr, and clusters formed 
during the last 500 Myr are excluded. In the right-hand panel, the dotted 
line gives a power law fit, which becomes shallower as the population ages. 



2011b). which is also shown by the assembly of the stellar halo 
between the second and sixth images of Fig. [7] However, this is not 
sufficient to prevent a radial trend of the peak mass, because dur- 
ing and after the final coalescence the disruption rate in the galaxy 
centre is still higher than in the outskirts. We have fitted log-normal 
functions to the cluster mass distribution as a function of galac- 
tocentric radius to quantify its radial variation. Together with the 
dispersion of the cluster mass distribution, this is shown in Fig. [6] 
for simulation 1ml 1 at two different times after the completion of 
the merger. The peak of the mass distribution of the surviving clus- 
ters in our simulations is initially (t = 3.8 Gyr) constant for radii 
R gc < 20 kpc, but at larger radii it depends on the galactocentric 
radius as M pea k oc i? g 7°' 8 , signifying a higher disruption rate in the 
galaxy centre. By t = 4.9 Gyr, the destruction of clusters on radially 
anisotropic orbits (see Fig. [3} has led to a shallower dependence of 
Mpeak oc RZ . The dispersion of the mass distribution follows a 
similar evolution. It is almost fully insensitive to R gc at t = 4.9 Gyr 
and steadily decreases with time. 

While the radial variation of M pea k evolves towards the radi- 
ally independent form that characterises observed globular cluster 
systems, there is no guarantee that this will still be the case at the 
present day. If the further evolution of the merger remnant is quies- 
cent, a radial trend of the characteristic mass might be reintroduced 
after a few Gyr, since the disruption of clusters proceeds more 
rapidly near the galaxy centre than at large radii. It may be pos- 
sible to erase this radial dependence again later on, for ins tance due 
to per turbations of the cluster orbits by minor mergers (cf. lOu et al.l 
l20ll . or any other perturbations that ma ke the galaxy potenti al de- 
viate from spherical symmetry (see also lFall & Z hang 200 1J) . The 
details of the further evolution of the cluster mass distribution and 
its spatial variation will depend on the cosmic environment, and 
cannot be followed in a major merger simulation. Whichever cos- 
mic conditions may govern the further ev olution of the clust er pop- 
ulation, the relative universality (see e.g. ljordan et alj |2007) of the 
globular cluster luminosity function suggests that these conditions 
should be commonplace if the cluster populations of our merger 
remnants are indeed the progenitors of present-day globular cluster 
systems. This is discussed further in Sect. [4] 



Contrary to the survival fraction of star clusters (see Fig.|4]and 
Eq.llOt. the peak mass and dispersion of the cluster mass distribu- 
tion does not exhibit any variation with starburst intensity, which 
is shown in the left-hand panels of Fig. [7J This contrast with the 
survival fraction arises because f surv is set at a different time dur- 
ing the merger than the shape of the cluster mass distribution. The 
cluster destruction rate is highest when the gas density peaks, at the 
height of the interaction and starbursts. The galactic environment at 
that stage is in fact so disruptive that star clusters can be destroyed 
irrespective of their masses. The slope of the cluster mass distribu- 
tion is therefore fairly constant at the times when the majority of all 
clusters is being destroyed. Only after the onset of galaxy coales- 
cence the frequency and strength of tidal shocks drop to a level at 
which the most massive clusters are able to survive, implying that 
low-mass clusters are preferentially disrupted. This leads to a flat- 
tening of the cluster mass distribution at the low-mass end. Figure[5] 
nicely illustrates this, as the slope of the mass distribution mostly 
changes after 1.8 Gyr. The shape of the mass distribution is thus 
set during the aftermath of the merger and therefore does not cor- 
relate with the cluster survival fraction, which is determined by the 
starburst intensity at earlier times. 

The right-hand panel in Fig. [JJ does show that the dispersion 
of the cluster mass distribution correlates with the peak mass as 
<7i gM oc M^ 33 . This relation arises because the high-mass end 
of the mass distribution does not vary much across all simulations, 
implying that low-mass cluster disruption sets both the peak mass 
and the width of the distribution. Present-day globular cluster sys- 
tems have a peak mass of abou t log (M pea k/Mp)) = 5-5.3 and a 
dispersion of ai og M = 0.4-0.5 jjordan et alj|2007l) . Extrapolation 
of the power law fit in Fig.|7Jto such peak masses yields too small 
dispersions (~ 0.2), indicating that the relation between dispersion 
and peak mass should flatten as cluster disruption proceeds. There 
are hints of a dependence of aw, M on Mpeak for old globula r cluster 
populations in the Virgo Cluster (e.g. Ijordan et al.ll2007t) . but in- 
deed the correlation is not as strong as we find for the cluster pop- 
ulations of young merger remnants. A direct comparison between 
the dispersion in our simulations and in the observed systems is ob- 
structed by the variation of the truncation mass among galaxies 
in the Virgo Cluster, which is not included in our models. 

3.5 Sensitivity of results to model assumptions 

For the simulations that are used in this paper, we have adopted 
certain initial conditions and made a number of assumptions. In 
this section, the presented results are tested for any dependences on 
two key assumptions. We verify the influence of (1) the evolution 
of cluster radii and (2) the adopted particle resolution. The former 
is important because the dominant cluster destruction mechanism 
(tidal shock disruption) scales with cluster density, while the latter 
should be checked to see if our simulations indeed resolve the gas 
structure that disrupts clusters (see the discussion in Sect. 12. 3t . 

Figure [8] shows the effect of the numerical resolution and the 
star cluster mass-radius relation on the SFH and the time evolu- 
tion of the total number of clusters over the course of our merger 
simulation lml 1. The numerical resolution is changed by a fac- 
tor of two up and down, while the added mass-radius relations im- 
ply a constant radius rj, = 3.75 pc and constant density rh/pc = 
3.75(M/10 4 M©) 1/3 . The SFH and the number of clusters are both 
normalised to their disc values to enable a straightforward compar- 
ison. The figure demonstrates that the relative change of the SFR 
and number of star clusters in the galaxy merger simulations with 
respect to isolated discs is not much affected by numerical resolu- 
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Figure 8. Influence of numerical resolution (left-hand panels) and the star 
cluster mass-radius relation (right-hand panels) on the SFR (top panels) 
and the number of clusters (bottom panels) in merger simulation lmll. The 
SFR and number of clusters are normalised to twice the value they have 
in each corresponding isolated disc simulation. The red solid lines show 
the benchmark result from Fig. [2] with a mass-radius relation and num- 
ber of particles as described in Sects. I2~2l and l2.3l In the left-hand panels, 
the dotted blue lines indicate the same simulation as lmll but with dou- 
ble the particle resolution, whereas the dashed blue lines denote the same 
simulation with half the particle resolution. In the right-hand panels, the 
dotted green lines represent the same simulation as lmll but using a con- 
stant cluster radius rj, = 3.75 pc independent of cluster mass, while the 
dashed green lines mark the same simulation with a constant cluster density 
r h /pc = 3.75(M/10 4 Mq) 1 / 3 . 



tion and the cluster mass-radius relation. At the end of the simula- 
tions, all differences are of the order of the statistical scatter. The 
only significant deviation seems to occur between t = 2 Gyr and 
t = 3.5 Gyr, when the number of clusters in the reference simula- 
tion is relatively unstable with respect to the other simulations. This 
can be traced to the stochastic variation of the SFR in the top panels 
(and thus the cluster formation rate). The top-right panel illustrates 
the statistical spread of the SFH over different realisations of the 
model, because the cluster mass-radius relation does not influence 
the SFH. In other words, the three shown SFHs are the result of 
identical boundary conditions. It implies that the dip in the SFH of 
the reference simulation at t ~ 2.4 Gyr and the subsequent rise of 
the SFR are stochastic. Due to the ongoing disruption of star clus- 
ters, this statistical variat ion is magnified in thei r number evolution 
(bottom panels, also see lKruijssen et al1l2011bh . We can therefore 
conclude from Fig.[8]that the time evolution of the number of clus- 
ters in galaxy mergers relative to isolated discs is not influenced by 
numerical resolution and the cluster mass-radius relation. 

The absolute number of clusters in the simulations does 
change for different mass-radius relations, because star cluster dis- 
ruption in galaxy discs and galaxy mergers is dominated by tidal 
shocks, for which the cluster disruption time-scale depends on the 
cluster density. However, when comparing the merger simulations 
to the corresponding isolated discs, this difference is offset by a 
similar change of the absolute number of clusters in both cases. The 
sa me holds for any possib le variation due to numerical resolution 
(cf. iKruiissen et alj|201 lbl Fig. 4). This validates the results pre- 
sented in Sects. |3~Tl43. 31 where we looked at the impact of galaxy 
mergers relative to isolated dies. 

We did not yet verify the impact of the cluster mass-radius re- 
lation on the cluster mass distribution (cf. Sect. 13.4b . The low-mass 
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Figure 9. The star cluster mass distribution at t = 4.8 Gyr for varying nu- 
merical resolution and mass-radius relations. As in Fig. [8] the red solid line 
indicates the reference simulation lmll, with blue lines denoting double 
(dotted) and half (dashed) the numerical resolution, and green lines repre- 
senting constant cluster radii (dotted) and densities (dashed). To account for 
the statistical scatter in Fig. [8] all distributions are normalised to the number 
of clusters in simulation lmll and clusters formed during the last 500 Myr 
are excluded. As in Fig. [5] the slope of the initial mass distribution is shown 
as a dashed line. 

(disruption-dominated) end of the cluster mass distribution attains 
a slope equal to the mass dependence of the disruption time-scale 
dFall & Zhan3l200ll) . This slope and the peak mass should there- 
fore be sensitive to the mass-radius relation (see the expression 
for the tidal shock disruption time-scale in Eq. [7}. In Fig. [9] we 
show the cluster mass function in the merger remnant of simula- 
tion lmll at t = 4.8 Gyr, for the same set of simulations as in 
Fig. [8] Again, the numerical resolution of the simulations does not 
strongly influence the result. It only affects the high-mass trunca- 
tion of the mass distribution because the particle mass limits the 
cluster mass (see Sect. l2.lt . However (and as expected), the clus- 
ter mass distribution does vary for different mass-radius relations. 
The question thus arises which range is possibly covered by the 
normalisation and exponent of the actual mass-rad ius relation. We 
include d an extensive discussion of this topic in IKruiissen et al.l 
J20 1 lbf) . in which we motivated our choice of the mass-radius rela- 
tion b y comparison to JV-body simulations of dissolv ing star clus- 
ters dBaumgardt & Makinol2 003; Kiip per et al.l2008l) . The adopted 
mass-radius relation (fl, = ih,4[M/10 4 Mq]*, with fi,,4 = 4.35 pc 
a nd S = 0.225) is con sistent with the mass-loss dominated regime 
of lGieles et alJJioTll) . However, the evolution of star clusters is ini- 
tially dominated by expansion and the mass-radius relation can not 
be expressed by a single power law for the entire cluster history. 
iGieles et al.l J201 lb show that the exponent of the mass-radius rela- 
tion varies from 5 = —0.25 during the expansion-dominated phase 
to S = 0.17 in the mass-loss dominated phase. This is close to our 
adopted value, which should thus be taken as an upper limit to the 
exponent throughout the evolutionary histories of star clusters. This 
rules out the green dashed line in Fig. [9] which shows the cluster 
mass distribution for 5=1/3. 

Before estimating lower and upper limits to the peak mass 
of the cluster mass distribution, the choice of the normalisation of 
the mass-radius relation th,4 should also be evaluated because it af- 
fects the peak mass. The typical radiu s of young star clusters of 
10 4 M Q is r h ~ 3.75 pc dLarsenll2004l) , which we adopted as the 
normalisation for the constant-radius and constant-density relations 
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in Figs. [8] and [9] whereas clusters in our reference simulation have 
fh,4 = 4.35 pc to match N-body simulations. The mean and median 
half-mass radiO of Galactic globula r clusters fall in the same range, 
with 4.3 pc and 3.0 pc, respectively (Harris 1996, updated 2010 ver- 
sion), showing little variation with cluster mass. Such little change 
after a Hubble time of evolution with respect to young clusters sug- 
gests that the adopted normalisations in Fig. [9] are representative 
of typical star cluster sizes for populations of any age. Combining 
this with the exclusion of the constant-density (S = 1 /3) simulation 
above, we can conclude that our reference simulation (5 = 0.225) 
provides a lower limit for the peak mass of 10 25 Mq. The green 
dotted line in Fig. [9] shows the mass distribution for 5 = 0, w hile 
observations of young clusters suggest S = 0.1 jLarsenll2004l) . A 
peak mass of 10 4 M can thus be interpreted as an upper limit 
to the peak mass in a merger remnant. Our conclusions should be 
specified further in a future work, by sampling the initial cluster 
radii from a certain (possibly cluster mass dependent) distribution 
function and including a more physically motivated description for 
the further radius evolution. 



4 CONCLUSIONS 

We have performed a numerical study of major mergers of 
comparable-mass disc galaxies, complemented with a sub-grid 
model for the ongoing formation and evolution of their star cluster 
populations. The simulations have been used to address the relative 
contributions of cluster formation and disruption over the course of 
a galaxy merger, and to investigate the potential formation of the 
metal-rich part of a globular cluster system. The main results from 
our model are as follows. 

(i) During a galaxy merger, the total number of star clusters de- 
creases. The increase of the star formation rate during merger- 
induced starbursts is compensated by a stronger increase of the 
cluster disruption due to tidal shock heating by dense gas. 

(ii) Although during certain episodes the destruction rate is high 
enough to disrupt clusters independently of their mass, over the en- 
tire course of a merger low-mass clusters are most strongly affected 
by the destruction. When considering increasingly massive clus- 
ters, their number decreases by a smaller amount during a merger. 
If the cluster sample is limited to massive and young clusters to 
mimic observational selection effects, the net destruction cannot be 
detected and changes to a transient increase of the number of clus- 
ters during the starbursts, in agreement with observational results. 

(iii) The relative decrease of the number of clusters is stronger for 
higher peak star formation rates, because the enhanced formation 
and destruction of clusters are both caused by the high gas density. 
This trend is weaker for higher masses and may be reversed above 
M ~ 1-3 x 10 5 M©, where a stronger starburst may produce more 
clusters than a weak starburst. In Eq. [10] we provide a generalised 
expression for the survival fraction of clusters as a function of the 
gas depletion time-scale, which reflects the intensity of the star- 
burst. 

(iv) The peaks in the cluster age distribution and star formation 
history can be offset with respect to each other due to the elevated 
cluster disruption rate at the height of a starburst. This offset can 



be as large as 200 Myr jKruiissen et alj|201 lbl) . which implies that 
while the cluster age distribution can be used to reveal the occur- 
rence of a starburst, it cannot necessarily be used to determine its 
time or duration. 

(v) The orbital kinematics of the star clusters in a merger remnant 
are isotropic within galactocentric radii of ~ 50-60 kpc due to the 
destruction of clusters on highly eccentric orbits. This value is sim- 
ilar to the result for the accretion of globular clusters from satellite 
dwarf galaxies JPrieto & Gn edin 2008), which shows that it may 
not be possible to distinguish between in-situ and ex-situ cluster 
formation based on solely the orbital (an)isotropy of the cluster 
population. 

(vi) The preferential destruction of low-mass clusters causes the 
power law initial cluster mass function to develop a peak at a mass 
of about 10 2 ' 5 Mq during the final coalescence of the galaxies. This 
is a lower limit, as the precise value depends on the relation be- 
tween cluster mass and radius, with the post-merger peak mass po- 
tentially reaching up to 10 4 Mq if the cluster radii are completely 
unrelated to their masses. The peak mass only weakly correlates 
with galactocentric radius due to the destruction of clusters on ra- 
dially anisotropic orbits, and (for the adopted mass-radius relation) 
increases by about 0.3-0.4 dex per Gyr after the completion of a 
merger. Young to intermediate-age (~ 2 Gyr old) merger remnants 
should display a peak in the star cluster mass distribution at about 
10 3 Mq due to the destruction of low-mass clusters (see Figs. [5] 
and[9}. 

(vii) After a merger is completed, the star cluster population is sim- 
ilar to what a young globular cluster system would look like. Firstly, 
the ejection of clusters from star-forming regions into the stellar 
halo produces a spatial distribution that is comparable to that of 
globular clusters. Secondly, the peaked cluster mass distribution 
is intermediate to that of young massive clusters and old globu- 
lar clusters. Thirdly, the high star formation rate during a merger 
is capable of producing clusters that are massive enough to survive 
for a Hubble time. 

Interestingly, the high disruption rate after (globular) cluster 
formation could lead to a mass distribution with a peak mass of 
10 3 -10 4 Mq on such a short time-scale that only little further dis- 
ruption is required to obtain the current peak mass of the globular 
cluster mass distribution. This would imply that even the subset 
of clusters on the widest orbits around their host galaxies would 
be able to reach it before the present day. If the ICMF of globu- 
lar clusters had a Schechte r-type truncation at the high-mass end 
dKruiissen & Cooperll20T l|), any further disruption would not yield 
an additional increase of the peak mas s because it t hen saturates 
at about 10% of the truncation masfl (Gieles 2009). This would 
thus lead to a 'universal' globular cluster mass function, indepen- 
dent of galactocentric radius and current galactic environment. The 
current peak mass of globular cluster systems through out the uni- 
verse indeed happens to be ~ 2 x 10 5 Mq (e.g. Ijordan et al.l 
120071) . arou nd 10% lower than the estimated truncation mass of 
their ICMF dKruiissen & Portegies Zwartll 20091) . 

The increased cluster disruption rate in galaxy mergers is 
driven by the high gas densities that also cause the burst of star 
formation. This indicates that the mechanism of enhanced disrup- 
tion is not necessarily constrained to major mergers, and can be 



These are determined by assuming that light traces mass, i.e. that the 
clusters are not mass-segregated. This may underestimat e the radii by up to 
a factor of 1.5, potentially depending on cluster mass (see lGieles et a l. 201 1 
for a discussion). 



9 This percentage applies if the ICMF of globular clusters followed a power 
law with index —2 below the trunc ation, and the mass dependence of the 
disrup tion time-scale is 7 ~ 0.7 jGielesI l2009l ; [Kruiissen & Po rtegies Zwart 
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generalised to any environment with a high gas density and a cor- 
respondingly high SFR. While major mergers may provide an effi- 
cient formation channel for globular cluster populations, they are 
not a prerequisite. Any extremely high-density, dynamically ac- 
tive, star-forming environment - be it in a starburst dwarf galaxy, 
during bulge assembly, in an unstable high-redshift disc or in a 
major merger - would cause the enhanced disruption of clusters 
at young ages. The clusters that eventually survive are charac- 
terised by a more quiescent evolution due to cluster m i gration and 



tensed by a more quiescent evolution due to cluster migration ana 
natural selection l Kruijssen & Portegies Zwart 20091; lElmegreenl 



120101: lElmegreen & Huntenl2010l : iKruiissen et alj|201 lbl) . Indeed 
the wide variety of galaxy types with remarkably similar globu- 
lar cluster mass distributions is hard to explain if cluster disruption 
is governed by the present-day environment, and suggests that the 
bulk of the disruption occurred at the epoch of globular cluster for- 
mation, when the host galaxies were likely more similar. A gener- 
alisation to all dense environments is supported by dwarf galaxies 
like Fornax, which has not undergone a major merger and yet har - 
bours a handful of globular clusters t Sha plevll 19391 ; iHodgdl 196lh 
that presumably formed in a starburst during the early formation of 
the galaxy. If such a generalisation to all dense environments in- 
deed holds, it would suggest that globular cluster populations may 
be the inevitable outcomes of the large starbursts occurring in the 
early universe. 
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